Sequence selectivity of the cleavage sites induced by topoisomerase I inhibitors: a molecular dynamics study

Topoisomerase IB (Top1) inhibitors, such as camptothecin (CPT), stabilize the Top1-DNA cleavage complex in a DNA sequence-dependent manner. The sequence selectivity of Top1 inhibitors is important for targeting specific genomic sequences of therapeutic value. However, the molecular mechanisms underlying this selectivity remain largely unknown. We performed molecular dynamics simulations to delineate structural, dynamic and energetic features that contribute to the differential sequence selectivity of the Top1 inhibitors. We found the sequence selectivity of CPT to be highly correlated with the drug binding energies, dynamic and structural properties of the linker domain. Chemical insights, gained by per-residue binding energy analysis revealed that the non-polar interaction between CPT and nucleotide at the +1 position of the cleavage site was the major (favorable) contributor to the total binding energy. Mechanistic insights gained by a potential of mean force analysis implicated that the drug dissociation step was associated with the sequence selectivity. Pharmaceutical insights gained by our molecular dynamics analyses explained why LMP-776, an indenoisoquinoline derivative under clinical development at the National Institutes of Health, displays different sequence selectivity when compared with camptothecin and its clinical derivatives.